Hard diffraction and the nature of the Pomeron 

J. Lamouroux*, R. Peschanski^^ , C. Royon-'- and L. SchoefFel^ 

We ask the question whether the quark and gluon distributions in the Pomeron obtained from 
QCD fits to hard diffraction processes at HERA can be dynamically generated from a state made 
of valence-like gluons and sea quarks as input. By a method combining backward Q^-evolution for 
data exploration and forward Q^-evolution for a best fit determination, we find that the diffractive 
structure functions published by the HI collaboration at HERA can be described by a simple valence- 
like input at an initial scale of order ~ 2.3 — 2.7 GeV^ . The parton number sum rules at the initial 
scale for the HI fit gives 2.1 ± .1 ± .1 and .13 ± .01 ± .02 for gluon and sea quarks respectively, 
corresponding to an initial Pomeron state made of (almost) only two gluons. It has flat gluon density 
leading to a plausible interpretation in terms of a gluonium state. 

1. Introduction 

Since years, the Pomeron remains a subject of many interrogations. Indeed, defined as the virtual colourless carrier 
of strong interactions, the nature of the Pomeron is still a real challenge. While in the perturbative regime of QCD 
it can be defined as a compound system of two Reggeized gluons [1] in the approximation of resumming the leading 
logs in energy, its non-perturbative structure is basically unknown. 

However, in the recent years, an interesting experimental investigation on "hard" diffractive processes led to a new 
insight into Pomeron problems. At the HERA accelerator, it has been discovered that a non negligeable amount of 
7*— proton deep inelastic events can be produced with no visible breaking of the incident proton. There are various 
phenomenological interpretations of this phenomenon, but a very appealing one (which indeed constituted a prediction 
[2]) relies upon a partonic interpretation of the structure of the Pomeron. In fact, it is possible to nicely describe 
the two sets of cross-section data from HI [3] and ZEUS [4] (after taking into account also a Reggeon component) 
by a QCD DGLAP evolution of parton distributions in the Pomeron combined with Pomeron flux factors describing 
phenomenologically the probability of finding a Pomeron state in the proton. Sets of quark and gluon distributions 
in the Pomeron following LO or NLO Q^-evolution equations are obtained [5], which successfully describe HI and 
ZEUS data sets separately^. 

The idea carried on in the present work is to find whether the Pomeron structure functions can be obtained from a 
standard DGLAP [6] evolution equation initialized by a valence-like input at some scale /x. By valence-like we mean 
an input distribution at low scale /x^ for which both the density of gluons and sea quarks remains finite when /? — > 0, 
where (3 is the energy-momentum shared between the constituants. If this is achieved, the number of gluons and 
sea quarks in the Pomeron is well-defined and finite. This may give an information on the non-perturbative origin 
of the Pomeron which can then be interpreted as a state made of constituent gluons and sea quarks with a given (3 
distribution (which reflects the energy-momentum sharing between the constituants and thus their interaction) and 
a given transverse size given by the scale /i. 

Our approach is inspired by the well-known GRV approach for the proton [7] (as well as its extensions to pion, 
photon) where the gluon, quark and antiquark distribution are obtained via dynamical parton generation from LO 
or NLO DGLAP evolution starting at low scale. In the case of the proton, taken as an example, the idea was in 
principle to generate all proton distributions from QCD radiation starting only with the three valence quarks. After 
a series of refinements due to precise data analysis [7], it happens necessary to introduce also valence-like antiquarks 
and gluons. In any case, the overall picture leads to the useful GRV set of structure functions widely used in QCD 
phenomenology. 

The main results of our analysis arc herafter summarized: 

i) The HI set of data is fully compatible with an initial valence-like set of sea {S = q-\-q) and gluon (G) distributions 
at a scale fij^^o 2.7 GeV^ (/i|o - 2.3 GeV^). 
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'^Both sets are compatible within the errors but a difference was observed between the two sets of structure functions corre- 
sponding to different trends in the Q^-dependence of the data [5]. This motivates the separation between the two sets of data 
for the analysis. 
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ii) Simple input valence-like distributions 

/35(/3,Q2 = M') = a,/3(l-/3) 

pG{p,Q^ = l,') = ag(3{l-p)'^ , (1) 

give rise to a good fit of data. The four parameters Ug^g, bg^g are given in table I with both statistical and systematic 
errors. 

ill). The values obtained for bg,bq are small, leading to flat initial parton density distributions for the Pomeron 

at the initial scale fi^. The exponential term in the usual paramctrizations of input quark and gluon densities [3,5], 
with an essential singularity at /3 = 1, is needed to compensate for the singularities of QCD matrix elements at low 
Q^. This explains the difference between the usual initial parametrizations for diffractive vs. total proton structure 
functions. 

iv) The ZEUS set of data, when parametrized using (1), leads to a worse fit, see table I. Even if these data have 
smaller weight in a global analysis, this difference deserves further study. 

v) Considering the obtained set of Pomeron structure functions, the input parton distributions of the Pomeron can 
be interpreted as those of a gluonium in a fundamental quantum state. Indeed, the gluon and sea quark number QCD 
sum rules (see formula (6) in section 4), correspond to a state of (almost only) two interacting gluons with (almost) 
flat parton density in /3. 



2. Backward evolution 



Our method can be decomposed in two steps. First, we consider the existing parametrizations of the Pomeron 
structure functions, which are not of GRV type. We perform a backward QCD evolution in order to see whether 
at the same small scale, the sea quark and gluon distributions can be compatible with the valence-like property. 
Then, if this investigation leads to a positive result, we perform a new QCD analysis of data, starting directly from a 
valence-like input. A consistency check is to verify that the new parametrizations of the Pomeron structure functions 
are compatible, within errors, with the initial ones. 

As an illustration of our method for generating dynamical parton distributions, let us first consider the backward 
evolution for the proton within the GRV scheme [7]. By construction, the parton distributions are obtained from 
valence-like quark, antiquark and gluon distributions at a small scale /x^. In the Mellin j-plane, the valence-like inputs 
correspond to Mellin transformed moments, or more generally continuous j-distributions, which remain finite when 
j — > 1. Indeed, the j = 1 moment, when it is well-defined and normalized to the j = 2 energy-momentum sum rule, 
defines the number of partons. It is constant for valence quarks for any Q^, while it is in general infinite for sea quarks 
and gluons, except eventually at the input scale /z^, if and only if the input distributions are valence-like. 

In the case of the total structure functions of the proton, this is what can be seen in Fig.l, where we have 
reproduced the NLO backward evolution back to /x^ of the valence, sea quark and gluon distributions from the GRV 
parametrization. It shows that all distributions which are infinite at Q > jj, become finite at j ^ 1 when Q ^ ji. 
In this instructive exercise we use the Mellin transform formalism for the (LO and NLO) Q^-evolution [8], which is 
particularly suitable for our purpose since DGLAP equations take a simple two-by-two matrix form in Mellin space 
and allow us to use an exact analytical solution for the Pomeron structure functions at all Q^. We use the same 
method for the diffractive structure functions. 

In our search for a valence-like input for the Pomeron, we now use the NLO backward evolution starting from the 
ansatze of Ref. [5]. Indeed, the key technical point of our analysis is the identification of an exact analytic Mellin 
transform of the input parametrizations commonly used [3,5] for the DGLAP evolution of the Pomeron structure 
functions. 

Starting from the parametrizations at = 3 GeV'^ [5], the sea quark distribution S'(/3, Q"^) = J2fiavorsil'^^)il^' 
and the gluon distribution G(/3, Q^) are parameterized in terms of coefficients C, 



(G,S) 



f3S{fi,Ql) 



^Cf)p, (2/3-1) 



i=i 



(2) 



where Pj{C) is the j member in a set of Chebyshev polynomials. Pi — 1, P2 — ( and Pj^i{() ~ 2(^Pj{{^) — Pj-i{Q. 
The parameters ' used in [5] are given in Table II. In formula (2), one takes the parameter a = .01 . 

Now, let the structure functions be expanded as J2t=o (^i^'^^ exp[a/ (/3 — 1)], with straightforward linear relations 
between d^'^'^^ and C^'^'^^; The Mellin transform reads [9]: 



2 



4 

F(^-^\j,Q') ^ J2 ^f''^ r(i+i-l) e-t Wi^_,-(a) , (3) 

where the Whittaker function Wi_^_^ i{a) can also be expressed in terms of the Meijer function G24 [9]. 

Once using expressions (3), together with the NLO or LO evohition scheme [8], it is straightforward to get the whole 
j-dependent parton distributions in any suitable range of Q^, cither for forward {Q > Qo) or backward {Q < Qo) 
evolution. 

For the Pomeron, looking for a valence-like input, we use the NLO backward evolution starting from the ansatze of 
Ref. [5] and look for the possibility of valence-like distributions in some range of ~ The analytic singularities of 
the Whittaker function in (3) are approximately cancelled by those of the evolution matrix elements. They arc both 
situated at j = 1 and the backward evolution induces a change of sign both for the sea and the glue distributions 
in the same range ~ /x^. 

The results (dashed lines) arc shown in Fig. 2 for the HI sea and gluon j-distributions and in Fig. 3 for ZEUS. 
Let us first comment the results obtained for the parametrizations of HI data. As shown in Fig. 2, the behaviour of 
the Mellin transformed distributions show a different trend for small values. The singularities present at j = 1 
in the initial parametrizations (3) interfere negatively with those present in the evolution matrix elements [8]. The 
key point is that they do so both for sea and gluon distributions in the same range of « 2 GeV^ below which a 
transition occurs. Aromid this value, the singularities in matrix elements overcome the initial ones. In the same figure, 
a qualitative exemple of forward evolution starting from a valence-like input, showing how the dependent matrix 
elements conspire to mimic the original parametrizations of [5]. This will be studied in detail and quantitatively 
confirmed in the next scc;tion. 

Note that the exponential term in formulae (2), which was necessary to describe the HI diffractive structure function 
data [3,5], has an essential singularity at /? = 1. It can now be understood as being induced by the compensation of 
the singularities of QCD matrix elements at low Q^. It is related to the flatness of the input parton distributions of 
the Pomeron^ that we find in our fits of the diffractive structure functions. 

Fig. 3 deals with the analysis of the backward evolution for the parametrizations of ZEUS data performed in paper 
[5] within the same scheme as for HI data. Quite interestingly, we did not find a range of below which both sea and 
gluon j-distributions meet a transition when j — > 1. In fact, while the glue distribution flips down at lower Q^, the sea 
keeps its singular trend. This result can be traced back to the difference seen in the /3-distributions of [5] . The gluon 
distribution is weaker than for HI, while the sea distribution remains larger at /3 ^ 0. In some sense, the backward 
evolution is more drastically driven by the sea than by the glue. This qualitative result of backward evolution will 
also be confirmed by the subsequent forward evolution analysis. 

3. Forward evolution 

It is well-known [3,5] that the diffractive structure function F^'-'^-', measured from DIS events with large rapidity 
gaps can be expressed as a sum of two factorised contributions corresponding to a Pomeron and secondary Reggeon 
trajectories. 

F^^^\q^(3,xv) = fv/^{xv)Ff{Q\p) + fu/^{xv)Ff{Q\l3) , (4) 

where xp is the fraction of energy of the proton flowing into the Pomeron or Reggeon. In this paramcterisation, 
can be interpreted as the structure function of the Pomeron and thus be expressed in a conventional way in terms of 
(3S{(}, Q^) and fiG{(5, Q"^). The same can be said for for F^, with the restriction that it takes into account various 
secondary Reggc contributions which can hardly be separated and whose structure functions are modelized from 
the pion ones. The Pomeron and Reggeon fluxes are assumed to follow a Regge behaviour with linear trajectories 
Q;p,R(i) = Q;p,]B,(0) + ccp pt, such that 

/p,R/p(a:p) = / 2a^,H(t)-l ^^ 

•'tcut 

where \tmin\ is the minimum kinematically allowed value of\t\ and tout = — 1 GeV^ is the limit of the measurement. 



^The situation is different for the total structure functions. The input distributions [7] are decreasing like powers when /3— >1. 



3 



The HI and ZEUS diffractivc data [3,4] were fitted using the simple formulae (1) for the sea quark and gluon 
distributions, with five free parameters: aq,bq, ag,bg and /x^, the initial scale of the forward evolution. In order to 
avoid a region where the diffractive longitudinal structure function is large, only data with y < 0.45 are included in the 
fit to find the Pomeron and Reggeon intercepts. Further more, for the Reggeized flux factors, we take the same values 
of the Pomeron and Reggeon trajectory intercept as in Ref. [5], namely ap(0) = 1.20 ± 0.09 and a_R(0) = 0.62 ± 0.03 
for HI, and ap{0) = 1.13 ± 0.04 for ZEUS (there is no need for Reggeons for ZEUS data). Only data points with 
> 3 GeV^, (3 < 0.65, Mx > 2 GeV, and y < 0.45 are included in the QCD fit to avoid large higher twist effects and 
the region that may be most strongly affected by a non-zero value of R, the ratio of the longitudinal to the tranverse 
diffractive structure functions. The QCD fit was performed both at leading order and next-to-leading order using the 
usual DGLAP [6] evolution equation. 

The fitted parameter values obtained for the HI and ZEUS data are given in Table I. The value obtained for the 
HI collaboration is quite good (x^=210.7 and x^ = 194.6 respectively at LO and NLO for 161 data points) and the fit 
result at NLO is shown in Fig. 4. We even note a good description of the data at high f3 which are not included in the 
fit (dashed lines in Fig. 4). 

In Fig. 5 and Fig. 6 are displayed the parton densities obtained with the LO and NLO fits. As expected, we find that 
the gluon densities is much larger than the quark one. The input distributions are characterized by a starting scale 
{fj,^ =2.7 GeV^) which is higher for diffractive structure functions than for the total ones [7] (.4 GeV^, see Fig.l). 

The error band corresponds to the systematic and statistical errors added in quadrature. Note that the errors 
are very small at low values of (3 by definition since we impose the behaviour at low (3 to be proportional to /3 in 
our parametrisation, see equation (1). As shown in the figures, the result is also compatible with the quark and 
gluon densities [5] found using the usual HI QCD fit using the input (2) with the parameters of table II. The scaling 
violations obtained using our parametrisation are given in Fig. 7. We note that our parametrisation leads to positive 
scaling violations for all /? values and flattens out at the highest values of /3, in good agreement with data. 

The results for the ZEUS collaboration are given in Fig. 8 for the NLO fit. The obtained (x^=44.6 and x^—b2.1 
respectively at LO and NLO for 30 data points) are two times worse than the one obtained by the usual QCD fit 
[5] . We note in Fig. 8 that our parametrisation cannot describe the data at low /3. The gluon are quark densities we 
obtain are respectively two times higher and smaller than the usual results of the fit to the ZEUS data. We are thus 
unable to reproduce ZEUS data using our parametrisation. 



4. Outlook: The Nature of the Pomeron 



As a brief summary of our study, the data on hard diffraction at HERA, as obtained from rapidity gap selection 
by the HI collaboration, are well described using a QCD evolution of parton distributions in the Pomeron, starting 
mainly from valence-like gluons (plus a small fraction of valence-like sea quarks) as an input at low scale. More details 
have been given at the end of section 1. In this last section we want to discuss the possible physical interpretation of 
this phenomenological result in terms of characteristic non-perturbative structures of the enigmatic Pomeron. Note 
that the similar question about the proton was the underlying physical motivation for the GRV parametrizations [7] . 

Using the HI fit results, and considering the j = 1 moments which are finite by definition of the model, it is possible 
to compute the number of quarks and gluons at the starting scale fj?: 

[dps 

■' = 0.13 ±0.01 ±0.02 



J(3d0{S + G) 

JdpG 
J (3 d(3{S + G) 



2.1 ±0.1 ±0.1. (6) 



We thus find around 2.1 for the gluon density and around 0.1 for the quark density, which is compatible with a picture 
of Pomeron made of two gluons at the initial scale of low . Note that this number has nothing of an initial input of 
our study, since only j = 2 moments have been constrained by energy momentum conservation to sum to 1 at all . 

The probability density of gluons G{(3, Q^) (and sea quarks S{(3, Q^)) gives other interesting complementary infor- 
mation. It is displayed in Fig. 9, for various values of Q^- It exemplifies the dynamical parton generation through QCD 
evolution with an increasing number of partons at small /? at the expense of those at high [3. In the top figure, we see 
the input probability density G(/3, jJ?) (and sea quarks S{(3, jj."^)) which is almost flat in f3. This striking feature means 
that the (almost) two gluons refferred to above, have an (almost) equal probability in terms of energy sharing. This 
feature indicates a state with two highly interacting gluons, which cannot be interpreted as independent constituent 
gluons. This is different from what has been observed for the proton, namely a constituent quark model for the proton 
or, also, from the peaked distribution [10] of quasi-free heavy quarks in mesons, for comparison. 

Summarizing the features of the obtained input distributions: 
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• The gluon distribution is largely dominant over the sea quark one. 

• The total j = 1 moments (see formula (6)) are near the value 2. 

• The probability densities are nearly flat in momentum fraction. 

These features are quite reminiscent of a gluonium, in a fundamental 0"*""*" or 2"'""'" states, where almost no qq excitations 
are present. Such a spectrum has been found^ in a 1 + 1-dimensional reduction of SU {N) gauge theories. 

The relatively high (compared to the proton case in GRV) value of the initial scale /x is also to be remarked. It is 
less stable than the other parameters, considering the variations between LO and NLO fits, but it stays in the range 
2.1 — 2.7 GeV^. This could be interpreted as an input Pomeron state having a rather large mass squared. Note the 
existence of 2++ gluonium states near-by in mass [11]. They also have been already discussed as possible candidates 
of gluonium states situated on the Pomeron trajectories [12]. 

As an outlook, It will be interesting to look for an analytic representation of Pomeron QCD structure functions for 
all in the spirit of the GRV parametrizations. On a more experimental ground, it will be very interesting to verify 
our conclusions with the data announced by the HI collaboration [13]. A first look seems encouraging and these data 
when publicly available will allow a better determination of the parameters of our input valence-like parametrization 
(1). We thus expect that they will give even more information on the elusive nature of the Pomeron. 
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TABLES 





dq 


K 




bo 




x' 


HI LO 
HI NLO 


0.13±0.01±0.02 
0.13±0.01±0.01 


0.40±0.02±0.13 
0.30±0.01±0.04 


1.9±0.1±0.1 
1.90±0.03±0.05 


0.20±0.01±0.02 
0.20±0.02±0.04 


2.3 ±0.08±0.50 
2.7±0.1±0.50 


210.7/161 
194.6/161 


ZEUS LO 
ZEUS NLO 


0.18±0.07±0.02 
0.18±0.02±0.02 


0.40±0.06±0.04 
0.30±0.02±0.07 


2.3±0.3±0.4 
2.30±0.02±0.30 


0.02±0.02±0.01 
0.01±0.01±0.01 


0.50±0.05±0.03 
0.90±0.02±0.05 


44.6/30 
52.1/30 



TABLE L Values of the parameters obtained at LO and NLO for the fit to HI data (two first fines), and ZEUS data (last 
two lines). The first error is statistical and the second systematic. 



parameters 


HI 


ZEUS 




0.18 ± 0.05 


0.41 ± 0.02 


dp 


0.07 ± 0.02 


-0.16 ± 0.03 




-0.13 ± 0.02 


-0.11 ± 0.02 


dp 

cP 
cP 


0.82 ± 0.40 

0.22 ± 0.06 
0.01 ± 0.04 


0.53 ± 0.30 

0.28 ± 0.25 
0.02 ± 0.11 



TABLE H. Parameters for quark and gluon input distributions in the Pomeron [5], see formula (2). The parameters are 
given at the initial scale Qo = 3 GeV^ . 
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FIGURES 




1 2 3 4 5 ^ 



FIG. 1. Backward -evolution of the proton j -distributions from GRV. The backward NLO evolution of the proton 
j-distributions calculated from the GRV parametrizations is shown at the initial scale I-l%lo ~ -4 GeV^ (dashed lines) and at 
= 12 GeV^ (continuous lines). Dark lines: Sea quarks; Dark grey lines: gluons; light grey lines: valence quarks. 



s 




FIG. 2. Q^-evolution of the Pomeron j-distributions for the HI case: comparison between backward and forward evolution. 
Dashed lines: Backward NLO evolution; Continuous lines: exemple of valence-like inputs. Curves from top to bottom are 
displayed for = 10^ 75, 12, 4.5 GeV^ and n%i^o- 
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FIG. 3. Backward Q^-evolution of the Pomeron j'-distributions: ZEUS case. The dashed hnes correspond to the Backward 
NLO evolution of the Pomeron j-distributions [5]. Same values as in Fig. 2. 
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FIG. 4. Result of the NLO QCD fit to the HI F^* data. The x2 is 194.6 /161 data points. Full Une: bins used in the QCD 
fit, dashed line: extrapolation to the bins not used in the fit. We note a very good description of the data over the full range 
using the simple parametrisation, see formula (1). 
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FIG. 5. LO parton strucure functions (grey bands) obtained with the fit to HI data. Left: quark density, right: gluon 
density, for four different values of Q^, namely 2.7 (the initial scale /x^), 4.5, 12. and 75 GeV'^. The error bands correspond to 
the systematic and statistical errors added in quadrature. The result is compared with the known [5] NLO DGLAP fit in full 
line for HI data, and dashed line for ZEUS data. 
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FIG. 6. NLO parton strucure functions (grey bands) obtained with the fit to HI data. Left: quark density, right: gluon 
density, for four different values of Q^, namely 2.7 (the initial scale 4.5, 12. and 75 GeV^. The error bands correspond to 
the systematic and statistical errors added in quadrature. The result is compared with the known [5] NLO DGLAP fit in full 
line for HI data, and dashed line for ZEUS data. We note a good agreement between our results and the usual HI gluon and 
quark densities. 
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FIG. 7. Scaling violations. The prediction for the structure function xp ■ F^^^^ {xw = 0.005) is presented as a function of 
in bins of (3, over the full range accessed. Data points are from Ref.[14]. 
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FIG. 8. Result of the NLO QCD fit to the ZEUS data. The x2 is 52.1 /30 data points. We note a bad description of 
the data especially at high /3. 
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FIG. 9. Sea quark and gluon probability densities. The distributions S{z, Q^) and G{z, Q^) are displayed for the NLO fit to 
HI data with the error bands as in Fig. 6. The LO fit (mean value) is also shown (continuous line). The top of the figure is for 
the valence-like input distributions at /i^. 
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